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Abstract 

Mixed action lattice calculations allow for an additive lattice spacing dependent mass renormal- 
ization of mesons composed of one sea and one valence quark, regardless of the type of fermion 
discretization methods used in the valence and sea sectors. The value of the mass renormalization 
depends upon the lattice actions used. This mixed meson mass shift is an important lattice artifact 
to determine for mixed action calculations; because it modifies the pion mass, it plays a central 
^ role in the low energy dynamics of all hadronic correlation functions. We determine the leading 

t— 5 order, O(a^), and next to leading order, ©(a^m^), additive mass shift of valence-sea mesons for 

\^ a mixed lattice action with domain-wall valence fermions and rooted staggered sea fermions, rele- 

vant to the majority of current large scale mixed action lattice efforts. We find that on the asqtad 
I— I improved coarse MILC lattices, this additive mass shift is well parameterized in lattice units by 

d A(am)^ = 0.034(2) — 0.06(2)(am7r)^, which in physical units, using a = 0.125 fm, corresponds to 

^ ^ A(m)^ = (291 zb 8 MeV)^ — 0.06(2)m^. In terms of the mixed action effective field theory param- 

os) eters, the corresponding mass shift is given by o^Amix = (316 ± 4 MeV)^ at leading order plus 

7^ , next-to-leading order corrections including the necessary chiral logarithms for this mixed action 

calculation, determined in this work. Within the precision of our calculation, one can not distin- 
guish between the full next-to-leading order effective field theory analysis of this additive mixed 
meson mass shift and the parameterization given above. 
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I. INTRODUCTION 



Recently mixed action (or hybrid) lattice QCD has become a commonly used method of 
computing hadronic observables with dynamical fermions in the chiral regime, as demon- 
strated by the rapid growth of mixed action calculations [HISlEllllElinilZlElEl [IDl [HI 
[T2I \T3\ [HI [151 [IB]- To compliment this progress, there have also been significant devel- 
opments in our theoretical understanding of mixed action (MA) lattice QCD through the 
use of MA effective field theory (EFT) which describes the low energy dynamics of the the- 
ory [I71[l8lll9l|2ni[2ll[22l[23l|2ll|25l|26ll^ Mixed action QCD allows the use 
of a numerically cheap variant of lattice fermions in the sea sector, where a significant part 
of the computational time is devoted, while using valence sector fermions that respect chiral 
symmetry at finite lattice spacing such as domain-wall [SU [321 [33] and overlap |3H [3S1 [3S] 
fermions. The use of these chiral fermions in the valence sector generally provides for both 
improved lattice spacing scaling and a simplification of the renormalization of lattice matrix 
elements [371 133 [321 UHl HI], in some cases eliminating potentially harmful power diver- 
gences. Additionally, the EFT extrapolation formulae for observable quantities computed 
with MA QCD have a form very similar to those for quantities computed with chiral fermions 
in the valence and sea sectors, in some cases an identical form through next-to-leading order 
(NLO) [271 [30], which have only small deviations from the continuum chiral perturbation 
theory formulae [121 [131 [B]. In other words mixed action lattice QCD provides most of the 
benefits of lattice fermions with chiral symmetry, while avoiding the large cost of using these 
fermions in the sea sector. 

The most popular mixed action scheme was developed by LHPC [H |S] and has been used 
in several recent calculations [Il[3l[ll[5ll6l[7l[Hl[ni HDl El HSl UHl [16] . In this approach domain- 
wall fermions are used in the valence sector while the fermionic determinant is represented 
by asqtad- improved [l5l US] fermions using the publicly available MILC configurations [17] . 
The MILC configurations are generated with rooted staggered sea fermions to recover the 
correct number of degrees of freedom in the continuum limit. There has recently been much 
discussion regarding the validity of the rooting procedure [IHl HH [SHI [SB [S21 1531 ISH [551 [5ni 
[521 |5H]| and a comprehensive review of the issues can be found in Ref . [52] • 

Independent of these issues related to rooted staggered fermions, mixed action calcula- 
tions come with their own set of unique challenges. At finite lattice spacing, MA theories 
violate unitarity. Also, hadrons composed of all valence quarks, all sea quarks, or a mix of 
valence and sea quarks each receive different lattice spacing dependent mass renormaliza- 
tions. However, these issues are addressed in a MAEFT framework [17] which is a natural 
extension of partially quenched chiral perturbation theory (PQxPT) [6OI [61]. The unitar- 
ity violations in a mixed action effective theory are identical in form to those in partially 
quenched effective theories which is transparent with the introduction of partial quenching 
parameters [231 [27] : these provide the means to quantify the unitarity violations at the EFT 
level. 

One particularly nice feature of MA theories is that the form of the leading order (LO) 
lattice spacing dependent additive mass renormalization for a mixed meson composed of one 
valence and one sea quark is universal, i.e. independent of the particular lattice actions used 
for the valence and sea sectors, and is described by one parameter, known in the literature 
as Cuix [20]. This parameter, although not physical, is an important lattice artifact to 
be determined for MA calculations; the long-range correlations of hadronic observables are 
modified by this parameter as the infrared behavior is dominated by the lightest particles 
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in the system, which include the mixed valence-sea mesons. There are some observables 
in which it has been shown this parameter does not appear even through NLO in the MA 
expansion [221 ESI , although this is generally not the case 



In this article we present our computation of the lattice spacing shift to the mixed 
pion mass, specific to domain-wall valence fermions and the rooted staggered coarse MILC 
fermions, relevant for the most commonly performed mixed action calculations at present. 
We perform this calculation by measuring the long-Euclidean time behavior of two-point 
correlation functions constructed with one domain-wall quark and one staggered quark, tied 
together to form mixed pions. As with all calculations using rooted staggered sea fermions, 
the validity, and perhaps usefulness of this work depends upon the the rooting procedure 
common to staggered lattice actions. By calculating the mixed pion mass, we are able to 
determine a numerical value of the MAEFT low energy constant, Cuix-^ We find that the 
LO description of the mixed meson mass shift in the MAEFT is insufficient to describe our 
computational results. Therefore, in Sec. [Tl| we determine NLO extrapolation formula for 



a mixed meson with chiral valence fermions and rooted staggered sea fermions. In Sec. Ill 
we present the details of our lattice calculation, highlighting some peculiarities which arise 
in the calculation of our mixed domain-wall-rooted-staggered meson correlation functions 



(the details are collected in the appendix). In Sec. IV, we analyze our results with the NLO 
extrapolation formula, finding good agreement. Additionally, we find a simple quadratic in 
parameterization of the mixed pion mass splitting provides an equally good description 
of this quantity, allowing for a user-friendly consumption of our results. We then conclude 
in Sec. El 



II. MIXED MESON MASSES TO NLO IN MAEFT 

We begin with a brief introduction of mixed action effective field theory to set our con- 
ventions and motivate our study of the additive mixed meson mass renormalization. The 
EFTs appropriate to describe the low-energy dynamics of mesons were first developed for 
Ginsparg- Wilson valence fermions and Wilson sea fermions in Ref. [T71 [TS] and for staggered 
sea fermions in Ref. [20]. As mentioned above, these theories are a natural extension of 
partially quenched chiral perturbation theory (PQxPT) [60l El] (with PQxPT recovered in 
the continuum limit, a — > 0, of MAEFT). The MAEFT is a dual expansion in the quark 
mass and the lattice spacing, as are all lattice EFTs, and is constructed in a two step pro- 
cess developed by Sharpe and Singleton [62]: one first constructs the Symanzik effective 
continuum theory for the corresponding lattice action [63l IM] and then one constructs the 
low energy effective theory corresponding to this Symanzik action. 

For the coarse MILC lattices, the LO lattice spacing effects, which begin at 0{g'^a?), 
are comparable to the LO quark mass effects, as seen by the splitting amongst the various 
staggered taste meson masses [65], which suggests a useful power counting is ^ 5^ where 

2 2 A 2 2 /1 \ 



^ We can unambiguously determine the LO mixed pion mass shift, o^Amix = a^lSCiviix//^, however the 
extracted value of Cmix will obviously depend on the choice of /, whether it is chosen at a given pion 
mass, or in the chiral but not continuum limit, or in both limits, etc. 
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At leading order in this power counting, the MA Lagrangian with GW valence quarks 
and any type of sea quarks is given by [HI [20] 

£ = :^str - ^str (m^S^ + SmJ) + (Usea - Kvs) , (2) 

where we use the normalization / ^ 132 MeV, and the sigma field contains the mesons, 

The matrices M and M contain bosonic mesons while x ^^nd contain fermionic mesons 
with one ghost quark or antiquark. For an explicit representation see for example Ref . |27] . 

In this article, we are interested in staggered sea fermions, for which the LO potential 
can be found in Refs. [66l l67] . The mixed meson potential has a universal form, regardless 
of the type of sea and valence fermions, and is given by 

W.s = CMixStr(r3ST3St) , (4) 

where the fiavor matrix T3 is a difference in projectors onto the valence and sea sectors of 
the theory, 

n = Vs~Vv = diag(-Iy , Is, -Iv) , (5) 

where ly and I5 are identity matrices which span the space of valence and sea fermions 
respectively. These mixed and sea potentials, Uys and Usea give rise to the LO additive 
lattice spacing dependent renormalization of the meson masses. The sea-potential splits the 
staggered meson spectrum into five different tastes at this order while the mixed-potential 
provides a flavor and taste independent mass shift to all valence-sea mesons composed of 
one valence quark, v, and one sea quark, s, giving rise to the LO parameterization of the 
spectrum 

^1,V2 = Boirriy, +"^^2) , 
fhl^ = Bo{my + rris) + a^Awix , 
^liS2,B = Boirris^ + m^J + a^A^ . (6) 

We are adopting the convention of using tildes to denote masses which receive additive 
lattice spacing dependent corrections |23]- The various taste splitting (a^A^) have been 
computed on the coarse (and flne) MILC lattices and are found to be 



a^As = , 

a^AA ^ (276 ± 1 MeV)^ 
o^At ~ (348 ±2 MeV)2, 
a^Ay ^ (404±2 MeV)2, 

a^A/ ^ (446 ±6 MeV)2. (7) 



The mixed meson mass correction. 



2 16 Civiix 



a AMix = a — — — , (8) 
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has not been previously determined, and it is this quantity in particular that we have 
computed.^ This is one of the more important unphysical parameters (lattice artifacts) in a 
given MA lattice simulation, because it modifies the mass of the mixed mesons, and therefore 
modifies the chiral behavior of all hadronic correlation functions, generally at NLO in the 
MAEFT extrapolation formulae [30] (for some observables this dependence does not appear 
until at least NNLO [231 ESI [27] ). 

To isolate this mixed meson mass renormalization we calculated the quantity 

Am^ = rhl^ - i(m^^ + m^^ ^) . (9) 

The method to calculate the mixed meson mass renormalization is not unique, however 
this particular choice has a few distinct advantages. First, as is obvious from Eq. ([6]), the 
LO quark mass dependence exactly cancels in this mass difference and thus the dominant 
contribution to this mass difference will be precisely the quantity of interest, o^Amix- Second, 
this choice is entirely unphysical not just at LO but to all orders in MAEFT: it must vanish 
in the continuum QCD limit. In particular, for degenerate valence and sea quark masses, 
this choice must be proportional to the square of the lattice spacing, and therefore is a clean 
and good way to directly determine the size of the MA lattice spacing artifacts. 



As we will show in Section IV, we find that there is non-trivial quark mass dependence in 
this mass splitting signaling a sensitivity to NLO effects. Therefore we not only determine 
the LO mass renormalization but also the NLO mass renormalization of the mixed valence 
sea mesons. This requires the NLO extrapolation formulae for all three meson masses, 
the valence- valence [20], the sea-sea [67] and the mixed meson. The NLO MA formula for 
the mixed valence- sea fermions has not been worked out before. It is a straightforward 
computation, similar to computing loop graphs in MAEFT as well as a staggered xPT [U5l 
EZl EHl EHl [70]. However, as this is not the focus of this work, we give here just one 
intermediate step and the answer. To be precise, we shall give the answer for a mixed meson 
mass composed of a domain-wall valence up-quark u, and a staggered sea down-qnaik I of 
arbitrary quark-taste t (not to be confused with meson-taste B). However, it will be clear 
how to generalize our expression for an arbitrary flavor non-singlet composition. First we 
provide an integral expression for the NLO loop correction to the mixed meson mass (as a 
Minkowski space momentum integral with regularization scheme R), 

t=l Fs=j,Lr B Fs=j,l,r 

+ A{k' - myg^^^^yik') + A{e - myg^^^^Jk')^ . (10) 

The first line of this expression contains all the physical contributions which survive the 
continuum PQ and QCD limits, while the other two lines contain only unphysical MA and 



^ In Rcf. '5^, this mass shift was estimated from a MA lattice calculation of the pion form factor using 
domain-wall valence and the coarse MILC sea quarks |3] by using the appropriate MA extrapolation 
formula. Unfortunately, only one lattice numerical result was in the chiral regime, and so only a crude 
estimate was made, o^Amix ~ (380 MeV)^. 
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staggering lattice artifacts. The flavor-diagonal momentum space propagators contain both 
"connected" (C) as well as "disconnected" (D) or haripin pieces, and standard expressions 
for these propagators can be found for example in Ref. ^T]. The meson-taste labels on the 
propagators, B = {5, A, T, V, 1} denote that the sea-sea meson masses in these expressions 
are of the taste B, with the mass renormalizations given in Eq. g. The expressions for the 
disconnected taste-l^ and taste-A propagators are slightly different from the PQ form and 
can be found in Ref. [67] . The explicit factors of ^ in front of the taste sums are necessary to 
account for the "rooting" of the staggered fermion determinant and can be understood either 
with the quark-flow picture or with the replica method [19| [73] . The summation over the 
four quark-tastes is denoted by a sum over t while the summation over the 16 meson-tastes 
is denoted as a sum over B, and the explicit flavor sums run over the sea-quark flavors, 
denoted j,l,r following the convention set in Ref. ^74j. Lastly, to determine the expression 
for a mixed meson of valence flavor V and sea flavor S, one simply makes the replacements 
in Eq. (10) u ^ V and / S, and then computes the given integrals in Eq. (10). 

We now provide the expression for the mixed pion mass to NLO in MAEFT of flavor 
ul. The NLO MA formula for the "charged" mixed pion mass in the degenerate light quark 
limit is given by (using dimensional regularization and the modifled minimal subtraction 
scheme conventional to xPT 



Bo{mu + mi) + a^AMix H 72 



+ Bo{mu + nil 



2L8(m) - L^ifi) 




vmip y 



In 



+ Go^At rhfp rp In 



m 



(11) 



In this expression, the mass of the sea-sea r^-octet fleld, which is of the taste-J variety, is 
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denoted m'j^^. The light quark partial quenching parameter is Af^ = mfij — m^. It turns 
out to be convenient to express the SU{3) breaking in the sea-quark masses in two ways, 

A^. = ^A^, = ^(<b--?.b)- (12) 

The new vector masses, m^^ and m^^ are the eigenvalues of a mass matrix composed of the 
quark mass contributions as well as the vector hairpin coupling a'^S'y contributions to these 
masses, for which the explicit form can be found in Eqs. (58) and (59) of Ref. [67], and 
similarly for the new axial meson masses and hairpin coupling as well. 

We now provide the NLO extrapolation formulae needed for our work with the mass 
tunings employed in most MA calculations to date [Il[3lll|5l|6l[3|8l[9l[T0l[IIl|12l|15l[T6],in 
which the domain-wall pion mass is tuned to the Goldstone staggered pion mass to within a 
few percent accuracy. This means we can treat all the quark masses and partial quenching 
parameters as equal. 
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a^A/, 



(13) 



The mixed meson mass renormalization defined in Eq. ^ is given through NLO by 

Am^ = (a'^Aui^ + a^L„4(/i)j + m^a^L„2„(/i) 



Mix 



(47r^)2 
^ 2a^AMix 



ml + a^Awix) In 



ml + a^AMix 



(m 



K 



a^AMix) In 



m]^ + a^A 



Mix 



(14) 



where x^^^ is a lengthy expression containing most of the chiral logarithms determined 



by Eqs. ([M]), Eq. (75) of Ref. [67| and for example Eq. (Bl) of Ref. [27]. We see 

from Eq. (14) that the mixed meson mass renormalization depends upon the LO mass 
renormalization a^AMix linearly, quadratically and also logarithmically. The other important 
thing to note is that x^^^ is a known function given the inputs of staggered lattice artifacts 
determined by MILC |;65j and that all the dependence upon the physical counterterms 
has exactly cancelled.^ Therefore, this function only depends upon three unknown terms, 
a^Anix, d'^La'i^fi) and a^La2m{fJ^)- With only one lattice spacing we can not independently 
determine Amix and a'^La* (/x) and so we treat these as one parameter in our fits which we 



report in Section |IV[ We now turn to the lattice calculations performed for this work before 
presenting the results of our analysis. 



^ This is only true for degenerate sea and valence quark masses, otherwise this mixed meson mass renor- 
malization depends upon one linear combination of Gasser-Leutwyler constants, the counterterms of 
XPT [m 113 E], Am2 cx i6(^°("'-'»^))' (ig _ 2Ls). 
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TABLE I: We lists the parameters of the MILC gauge configurations and domain-wah propagators 
we used in this work. With the mass mOlO lattices, we used propagators with both anti-periodic 
boundary conditions as well as Dirichlet boundary conditions. For this work, we generated an 
equal number of staggered propagators with masses equal to those in the configurations. 



Ensemble 
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arus 


2064f21b676m010m050 


0.010 


0.050 


2064f21b676m010m050 


0.010 


0.050 


2064f21b679m020m050 


0.020 


0.050 


2064f21b681m030m050 


0.030 


0.050 


2064f21b683m040m050 


0.040 


0.050 


2064f3b685m050 


0.050 


0.050 



dw f 

anil 


10^ X arrires"' 


# propagators 


0.0138 


1.552 ±0.027 


171 Full-V 


0.0138 


1.552 ±0.027 


447 DBC 


0.0313 


1.239 ±0.028 


222 DBC 


0.0478 


0.982 ± 0.030 


564 DBC 


0.0644 


0.834 ± 0.024 


349 DBC 


0.0810 


0.726 ± 0.030 


416 DBC 



"Determined by the LHP Collaboration. 

III. LATTICE CALCULATION 

A. Lattice actions and parameters 

Here we are computing the additive mixed meson mass renormalization relevant to the 
lattice actions used in a number of recent mixed action lattice computations [31 IH [5l El [H 
El El Uni [m [121 [16] • These calculations use domain- wall valence quarks with Nf = 2 + 1 
asqtad-improved [l5l |l6] quarks included in the publicly available MILC configurations |17] . 
In the valence sector, HYP-smearing [751 [ISl IZH] of the gauge field entering the domain- 
wall fermion Dirac operator has been used in order to reduce the explicit chiral symmetry 
breaking present at finite 5th dimensional extent. We worked at a lattice spacing a ~ 0.125 
fm keeping the strange quark mass fixed near its physical value, arUs = 0.05 and used a range 
of degenerate light quark masses. See Table [T] for details. The domain-wall propagators we 
used were generated by the LHP and NPLQCD Collaborations using Dirichlet boundary 
conditions (DBC) and reducing the time extent of the original MILC lattices from 64 to 32. 
To study the systematics introduced by Dirichlet boundary conditions we additionally used 
a set of full volume anti-periodic domain-wall propagators generated on 171 of the mOlO 
lattices. In all cases the domain-wall fermion action had an extra dimension of Lg = 16 and 
a domain- wall height of Mg = 1.7. The domain- wall fermion quark mass was tuned such 
that the domain-wall pion mass is equal to the staggered Goldstone pion mass within a few 
percent [HIS]. To construct the mixed mesons, we computed staggered propagators on these 
lattices with masses set equal to the corresponding sea quark masses. For computational 
simplicity, we used the same boundary conditions for the staggered fermions as was used 
on the domain-wall and fermions. All propagator calculations and meson interpolating 
operators were created with the Chroma software suite [79] . 

B. Meson interpolating operators and two-point correlation functions 

In order to compute the mixed meson mass renormalization defined in Eq. ([o]), we use 
a source for each of the domain-wall pion, the staggered Goldstone pion and the mixed 
domain-wall staggered pion. For the domain-wall pion we use the standard pion interpolating 
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operator 

7r(x) = ipdwix) 75 ipdw{x) . (15) 
To create the staggered pion and the mixed mesons, we use a "Wilsonized" staggered quark, 

ij:{x) = {n{x)cr x{x), (16) 

where a is a Dirac-spinor index, C can be chosen as a constant Dirac-spinor, x(x) is the 
single component staggered fermion field and 

4 

is the Kawamoto-Smit multiplicative phase factor [8QJ which accomplishes the Susskind spin- 
diagonahzation |811 |82l |83] relating staggered fermions to naive fermions, and we are using 
Euclidean Hermitian gamma matrices, 

{7;., 7.} = 25;... (18) 

This allowed us to use the existing Chroma software supplemented by a minimal amount of 
new software to create the staggered propagators as well as the staggered meson and mixed 
meson source and sink operators. For example, the staggered Goldstone pion can be created 
with the interpolating operator 

7r5(x) = tpsix) 75 (a;) 

oc e(x) xix) xix) , (19) 



which is known to have the quantum numbers of a spin zero pseudo-scalar meson [HI] and 
the phase is e(x) = (— l)^^^^^/'^). In this work, we not only look at the mixed pseudo-scalar 
correlator created with the interpolating operator, 

vr^, (x) = ipdu, (x) 75 i^s (x) , (20) 

and its Hermitian conjugate, but we also study the two-point correlation functions of all the 
mixed meson sources created with the interpolating operators 

Tilsix) = ipdw{x)Tipsix) , (21) 
with 

r = {1, 75, 7m, l5li,,7ti7u} , (22) 

and /i 7^ z/. The reasons for this will become clear below. This construction is very similar 
to that used for heavy mesons with staggered light quarks [HS]. There is a point we should 
note here. The virtual mixed mesons which propagate in internal loops of of mixed action 
hadronic correlation functions are composed of one domain-wall or overlap valence fermion, 
and one rooted staggered sea fermion. To determine the mixed meson mass, we construct a 
correlation function with a domain-wall fermion and a non-rooted staggered fermion. One 
might be concerned that the non-locality of the rooted staggered action [M] invalidates this 
method of determining the mixed meson mass. However, we need make no more assump- 
tions about the vanishing non-localities in the continuum limit than are necessary for the 
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construction and analysis of any correlation function with staggered valence fermions on 
staggered sea fermions. For example, this is how the MILC Collaboration determines the 
masses of the various taste mesons which propagate in virtual loops of correlation func- 
tions constructed with only the taste-5 mesons [55] • Furthermore, it is in this sense that a 
staggered lattice action in which the valence and sea quark masses are degenerate is still a 
partially quenched theory ^9] . 

The zero momentum, two point correlation function created by the domain-wall pion 
source is given by"^ 

a(t) = ^(7r(x,t) 7rt(0,0)) 

X 

= J2 ^ncosh (mn{t - , (23) 

n 

on the full periodic volume^ and on the lattices with Dirichlet boundary conditions with the 
replacement cosh — > exp. The staggered meson interpolating operator creates both states 
which are "straight" and those which oscillate in time, given by 

a,(t) = $^(vr5(x,t)4(0,0)) 

X 

= ^ A,cosh(^m„(t - + ^5„,(-l)"*cosh(^m„/(t - Nt/2)^ , (24) 

n n' 

with fit = t/a (and cosh — >• exp for DBC). For the staggered states created with the staggered 
meson interpolating field, the states which oscillate in time also have opposite intrinsic 
parity as compared to the straight states [HI], and are therefore significantly heavier than 
the lightest Goldstone pion. 

The mixed meson correlation functions also contain states which have a time-oscillatory 
behavior. The mixed meson correlation functions created with the arbitrary mixed meson 
interpolating operators, Eq. (|2T|, are given in the full volume by 



X 

oo 

= Y,[Al + (-l)"'i?g cosh(m„(t - iV,/2)) . (25) 

n=0 

In Appendix |X| we use the symmetries of the mixed domain-wall valence and staggered sea 
lattice action to show that the lightest mass state which dominates the long time behavior 
of these mixed meson correlation functions are always pions, such that for all F considered 
in Eq. (|22|), 

mo = rh^s ■ (26) 



^ Higher energy states in the sum will not necessarily correspond to single particle mass eigenstates but as 
we are only interested in the ground state we use m„ instead of En to represent the energy of a given 
mode in the sum. 

^ Anti-periodic boundary conditions for the quarks gives rise to periodic boundary conditions (PBC) for 
pions. 
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We find tliat the situation with Dirichlet boundary conditions in time is different. In this 
case we find that the long Euclidean time behavior of the mixed meson correlation function 
is given by 

Clfit) ^ <^e-'^-* + {-irsZ^e-^-' (27) 

with rhys 7^ M^s- For fermions with anti-periodic boundary conditions in time, there are 
two interpolating operators which couple to the mixed pions that are related by the time- 
doubling symmetry of the staggered action, in fact requiring m^s = M^g, as we observe with 
our full- volume calculations (see Appendix |A] for details). The Dirichlet boundary conditions 
break this time-doubling symmetry, which in position space, is the time-shift symmetry. We 
hypothesize this symmetry breaking is responsible for the observed lifting of the degeneracy 
between straight and oscillating states. In Appendix A4| we collect our results for rhys and 



Mys, demonstrating that both states must be pions, as they are too light to be anything else. 
Furthermore, one notices that the splitting in the square of these masses is independent of the 
pion mass, within errors, indicative that the splitting is due to a lattice spacing correction, 
consistent with expectations from MAEFT. As an important test of systematics, on the 
mi = 010, rris = 050 coarse ensembles, we have performed our calculations using a set of 
full-volume propagators satisfying anti-periodic boundary conditions in time, as well as the 
half volume calculation with a Dirichlet wall. In this case we find that 

m^^^ = 0.283(3) for time Dirichlet BCs , 

rhj,^^ = 0.285(3) for time anti-periodic BCs . (28) 

Given that the extracted mixed meson mass is the same within errors, the issue of the 
boundary conditions is tangential to the focus of this work, and we do not address it further 
in the main text. In Appendix |A 4 we provide further details of the connection between the 



breaking of the time-shift symmetry and the related change in staggered quark taste. 



IV. ANALYSIS AND RESULTS : MIXED MESON MASS RENORMALIZATION 

To determine the meson masses, we use a time-correlated x^-minimization fitting routine. 
We perform both standard and jackknife fitting methods. We fit the leading exponentials 
(cosh functions) for the domain-wall and staggered mesons to the two point correlation 
functions and the leading two exponentials for the mixed meson two-point correlation func- 
tions. For the full volume mixed meson correlation functions, we use both a single cosh fit 
with a straight as well as oscillating amplitude and we also perform the fit allowing for two 
non-degenerate masses. Explicitly, we fit the domain-wall pion correlation function with 

A^cosh(m^{t - Nt/2)^ for PBCs 

j^D^-m^t _ ^29) 

For the staggered Goldstone meson correlation function, we fit with an identical form, 

A^gcosh (m^^ {t - for PBCs 

^z)g-m^,i forDBCs. (30) 
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TABLE II: Pion masses and A(am)^ = (am^'*)^ — ^[(07717^)^ + (0771,^5)^] • We present our results 
for the masses determined with the fitting routine described in Sec. |IV[ We compute the additive 
mixed meson mass renormahzation both with the staggered meson masses we determined for this 



work as 


well as 


by inputting the sta^ 


;gered meson 


masses from Ref. |47j. 


see text for details. 


Ensemble 


pion 


result 


fit rang 


e 


/ d.o.f. 


Q 


m050 


DBC 




0.500(2) 


6-13 




4.9 / 4 


0.29 


m040 


DEC 




0.457(2) 


6-13 




6.4 / 4 


0.17 


m030 


DBC 


ml' 


0.412(2) 


9-13 




1.1/1 


0.29 


m020 


DBC 


ml' 


0.353(3) 


9-14 




0.54 / 2 


0.76 


mOlO 


DBC 


ml' 


0.283(3) 


5-12 




5.0 / 4 


0.29 


m050 


DBC 




0.477(1) 


6-15 




2.9 / 8 


0.94 


m040 


DBC 




0.433(2) 


5-11 




7.1 / 5 


0.22 


m030 


DBC 




0.375(1) 


7-14 




7.8 / 6 


0.25 


m020 


DBC 




0.315(1) 


6-14 




3.9 / 7 


0.79 


mOlO 


DBC 




0.224(1) 


D— iO 




D.i / O 


U.D4 


m050 


DBC 




n A^^(9\ 


1 in 




17/0 


n /I Q 
U.4o 


m040 


DBC 




n dd()(9\ 


5-8 




1.2 / 2 


0.55 


m030 


DBC 




0.381(2) 


6-8 




0.16 / 1 


0.69 


m020 


DBC 




0.319(2) 


4-6 




0.22 / 1 


0.64 


mOlO : 


Full-V 


ml' 


0.285(3) 


6-31 




20 / 23 


0.65 


mOlO : 


Full-V 




0.2235(05) 


7-31 




30 / 23 


0.14 


mOlO : 


Full-V 




0.2248(03) 


13-32 




17 / 18 


0.51 


Ensemble 


mOlO(Full-V) 


mOlO 


m020 


m030 


m040 


m050 


A(am)^ 


0.031(1) 




0.024(2) 


0.026(2) 


0.018(2) 


0.019(2) 


A(am) 


Itag. input 0.031(1) 


0.030(2) 


0.027(2) 


0.028(2) 


0.020(2) 


0.019(2) 



For the mixed meson correlation functions with DBC we fit the masses with the two expo- 
iiGnt 1 clIs 

Ale-^-^ + (-l)-5le-^^^^* (31) 
For the full volume lattices we fit to both 



as well as 



+ (-1)"'5.„J cosh(mr(t - N,/2)^ 
A.,,cosh(mr(t - iV,/2)) + (-l)"*5.„^cosh(M-(t - iV,/2)) 



(32) 
(33) 



to explicitly verify the degeneracy of the mixed pion masses, ml' and M!^' with the full 
volume propagators. The results of these fits are collected in Table [TTj The central values 
are determined by standard error analysis and the error bars by a jackknife analysis.^ 



^ In Appendix A 3 we present the results of fitting the various mixed meson two-point correlation functions 
and in Appendix | A 4| we present the results of the mass splitting between the two light mixed pions with 
DBCs. 
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We find that the Dirichlet boundary conditions cause significant difficulties in extracting 
parameters from the staggered pion two-point correlation functions/ As can be seen from 
Table [IT], we are only able to determine the staggered pion mass over a very short range 
of time as compared to the much greater time lengths used in typical staggered lattice 
computations with anti-periodic in time boundary conditions, for example in Refs. [HI 
The difficulty arises from a combination of factors. Because we are calculating pion two- 
point correlation functions, the backwards propagating states from the Dirichlet wall are 
light enough that they contribute to the signal to such an extent that they pollute the 
extraction of the forward propagating signal we are interested in even at short times from 
the source. This problem is then amplified by the time-oscillatory nature of the staggered 
correlation functions, becoming more problematic for lighter quark masses. We were not 
able to determine a mass of the staggered Goldstone pion (taste-5) on the mOlO lattices. 
On the m020-m050 lattices the staggered pion mass we determine is consistent with those 
in Ref. |17] within our errors. For both the mixed and domain-wall pion, we find that the 
masses we calculate on the DBG lattices are consistent with those determined on the full 
volume lattices, as is seen in Table |TT| 



A. Determination of mixed meson parameters 

To determine the mixed meson mass renormalization we perform a jackknife analysis of 
the difference 

A(am)2 = {aml'f - ]^{am^f - \{am^,f . (34) 

The results of this analysis are collected in Table |TTj For the reasons discussed in the previous 
section, when we calculate this quantity we use the values of the staggered meson determined 
in Ref. [Hj. In this case, the uncertainty in the input staggered pion masses is much less than 
the errors we determine for the mixed and domain-wall pions and therefore can be ignored. 
We determine this additive mixed meson mass renormalization using both the DBG as well 
as the full volume propagators on the mOlO lattices. As is clear in both cases, there is an 
upwards trend of the mass renormalization with lighter quark mass, a signal of NLO effects, 
(9(a^m^). To analyze these results we use several different fits. First we fit the NLO MA 



extrapolation formula given in Eq. ( 14 ) with a ^^-minimization routine to the lightest three 
data points, the lightest four and all five data points given in Table |TT| for both sets of mOlO 
results. In order to perform these fits, we use the values of the staggered taste splitting 
listed in Eq. ([T]) which are taken from Ref. [65]. We also input the values of the axial and 
vector hairpin couplings for which updated values can be found in Ref. [86j. These fits are 
all performed at the scale afi = 1 and with the choice af = 0.103. The difference between 
using this value for af, the value in the chiral limit, or the value at another pion mass is 
higher order in the MAEFT than we are working. Ghoosing a different renormalization scale 
will adjust the parameters of the MA fit but leave the mixed meson mass renormalization 



^ As a reminder, most of the domain-wall propagators we used came with Dirichlet boundary conditions. 
We did not have the computing time to make our own, and so we borrowed them from the LHP and 
NPLQCD collaborations. The choice of Dirichlet BCs for the the staggered propagators we did compute 
was one of convenience. Clearly, this would not be a good choice for performing staggered lattice QCD 
calculations. 
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TABLE III: Determination of mixed action and quadratic fit parameters. The upper table contains 
the fits in which we use the DBC propagators for the mOlO lattices. The lower table contains the 
fits in which we use the full volume propagators for the mOlO lattices. In both cases, we input the 
staggered pion mass values from Ref. |37|. 





Fits to NLO MA formula, Eq. (14) 


Fits to quadratic formula, Eq. ( 35 ) 


DBCs 


a2(a2AMix) a^{o?L^2^) Q 


a^(a^5Mix) a2(a2£^2^) Q 


m010-m030 
m010-m040 
m010-m050 


0.0383(15) -0.054(31) 0.30 
0.0396(11) -0.087(19) 0.25 
0.0394(10) -0.084(13) 0.42 


0.0309(29) -0.029(28) 0.33 
0.0336(22) -0.064(17) 0.19 
0.0336(19) -0.063(12) 0.35 


Full Volume 






m010-m030 
m010-m040 
m010-m050 


0.0394(10) -0.074(25) 0.19 
0.0401(08) -0.094(15) 0.25 
0.0399(07) -0.090(11) 0.40 


0.0332(22) -0.048(25) 0.22 
0.0350(17) -0.072(14) 0.20 
0.0347(15) -0.070(10) 0.36 



unchanged as this is a scale independent quantity. As an example, the NLO MA fit to the 
lightest four points calculated with the propagators satisfying DBCs at the scale a/i = 0.5 
gives a?{a?^Ma) = 0.0392(16) and a'^{a'^La2m) = —0.081(13) which directly verifies that the 
scale dependent contribution from the term a'^{a^La4.{fi)), in Eq. (14), is small compared to 



the LO term, a^AMix-^ This also demonstrates that the size of these unphysical counterterms 
are fairly insensitive to the renormalization scale. We present the results of these fits in 
Table III We additionally use a x^-minimization to fit the same sets of data to a quadratic 
form, 

A(am)^ = a^(a^(5Mix) + (fl"^7r)^o^C2m , (35) 



for which we also present the results of our analysis in Table III 

Fitting to these two different functions gives rise to different values of the unknown 
(counter)-terms, as is seen in Figure [l|^a), where we display the central values, 68% (dashed 
ellipses) and 95% (solid ellipses) confidence intervals for four different fits. The error ellipses 
on the left of Fig. ^a) correspond to fitting the data in Table [ll] to the quadratic form, 
Eq. (35). The larger ellipse (blue) is fit to the lightest four data points including the mOlO 
point calculated with the DBC propagators. The smaller ellipse (purple) corresponds to a 
fit of the lightest four points but using the results of the full volume propagators for the 
mOlO calculation. The right two error ellipses correspond to fitting our results of the lightest 



four mass points to the NLO mixed action formula, Eq. (14) for the results obtained with 
DBC propagators (larger red) and with the full volume mOlO propagators (smaller green) 
respectively. As is seen by these error ellipses, there is no significant difference in using 
either the results calculated with the DBC propagators or the full volume propagators. 

In Figure [l](6) we plot the resulting mixed meson mass renormalization, Eq (34) and cor- 
responding 68% confidence error bands resulting from the quadratic and NLO mixed action 
fits to the lightest four data points with all results calculated with the DBC propagators, as 
well as the results of our calculation listed in Table [11} For the mOlO point {am-,r — 0.2235), 
the lower (blue) point is calculated with the DBC propagators and the slightly higher (red) 



The a^Ajviix term is scale independent. 
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(a) (&) 

FIG. 1: (color online) This is a plot of the central result of our paper. In Figure [T]^a) we present 
the central values, 68% (dashed lines) and 95% (solid lines) confidence intervals for four different 
fits, see Sec. IV A for details. In Figure [T][| 6), we plot both the results of our computation as well 
as best fit curves and 68% confidence bands as a function of the domain- wall pion mass, am.,^. The 
darker curve and error band are the results of the quadratic fit to the lightest 4 data points and 
the lighter curve (and error band) represent the results of fitting to the NLO MA formula. The 
dashed curves represent the same fits respectively but to the lightest three data points. 



point is calculated with the full volume propagators. The darker (black) curve and error 
band correspond to the quadratic fit while the lighter (purple) curve and error band corre- 
spond to the NLO mixed action fit. We additionally plot the resulting dashed curves for 
the same quadratic (black dashed curve) and NLO mixed action (lighter dashed curve) fits 
to the lightest three data points. 

As can be seen from Figure [l](6), at the 68% confidence level, these fits are all consistent 
(which also includes all fits listed in Table III). Furthermore, for the lightest three data 
points they are equivalent, with very similar 68% error bands. In order to distinguish 
between these two curves one would need at least an order of magnitude increase in statistics. 
Therefore we conclude that within our precision, one can parameterize the mixed meson 
mass renormalization with either the NLO mixed action formula, Eq. (14), or the quadratic 
formula, Eq. (35): ^ 



A{amy = 0.040(1) - 0.09(2)(am^)^ + Xs 



NLO 



or 



A{amY = 0.034(2) - 0.06(2)(am^) 



(36) 
(37) 



where X2^^ is slightly different from x^^^ ^ind can be determined with Eq. (14). In these 



,NLO 



expressions, we have chosen the fits to the lightest four points for which all the mass renor- 



^ We note that in Ref. [21] an estimate of this quantity was made by comparing the MA extrapolation 
formula of the pion form-factor to a MA lattice simulation of this quantity, in which the mass renor- 
malization was estimated to be a^Ajviix ~ (380 MeV)^. In addition, an independent preliminary lattice 
calculation similar to this work is consistent with our results |87j . 
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malizations were computed with the DBC propagators and the staggered pion masses were 
taken from Ref. We re-emphasize that these are not distinguishable to the resulting 

functions in which the full volume propagators are used in the calculation. Using Eq. 



we can then determine Cmix from Eq. (36), however, we emphasize that what is more im- 



portant is the determination of the mass renormalization itself, Eq. (34). We find (recalling 

AMix ^ (708 MeV)^ , 

CmIx (or a TV/r \T\4: 



p 



(354 MeV)^ , (38) 



and choosing / = 132 MeV, 



Cmix ~ (255 MeV)^ . (39) 
It will be interesting to compare these same parameters determined on the fine MILC lattices. 



V. DISCUSSION 

We have calculated the additive mass renormalization of mixed mesons composed of one 
domain-wall valence fermion and one staggered sea fermion for the asqtad-improved [l5l 
l¥6j . publicly available MILC lattices [47|| which are relevant to a large number of recent 
mixed action lattice calculations P[3lll|5l|6l[71[8l[9l[l0l[IIl[T2l[T5l[l6]. This is an 
important unphysical parameter, or lattice artifact to determine for a given mixed action 
lattice computation. It modifies the mass of the mixed valence-sea pions and therefore plays 
a central role in the low energy dynamics of all hadronic correlation functions, entering 
generally all extrapolation formula at the one-loop level determined with the appropriate 
mixed action effective field theory [30]. 

The chiral extrapolations of the observables computed in this mixed action scheme [U 
[6] have not included this mass renormalization as it has not been previously determined. 



Now one can make use of either Eqs. (36) or (37) to determine the error in neglecting 
this mass renormalization in performing chiral extrapolations and hopefully improve the 
uncertainty in many of the observables calculated to date [HEllllEliniElIHlEl [IHl Ell [121 
[T^ [TB] . In particular it is known that the pseudo scalar decay constants depend upon the 
mixed meson masses at next-to leading order [20j. Additionally, it is known that in mixed 
action lattice calculations, the J = 3/2 Kit scattering length [27j, the pion form factor [2^ . 
baryon masses [22], nucleon axial charge [29], nucleon twist-2 matrix elements, the nucleon 
to delta electromagnetic transitions and the neutron electric dipole moment [30] to name a 
few, all depend upon the additive mixed meson mass renormalization as they depend upon 
the mixed meson masses at the one loop level in the mixed action effective field theory. 
Therefore, with this additive mixed meson mass renormalization known, one can perform 
mixed action chiral extrapolations of all these quantities, as well as any others computed with 
domain-wall valence fermions and the asqtad-improved staggered sea fermions at this lattice 



spacing, a ~ 0.125 fm. Clearly, Eq. (37) is the most useful for including this systematic 
effect in these extrapolations. To be used in the type of mixed action calculations this 
work is addressing [Il[3llll|5l[6l[7l[8l[9l[l0llIIl[I21ll5lll6], one takes a given mixed 
action extrapolation formulae [30], and wherever there is a mixed valence-sea meson in the 
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expression, one makes the replacement (in lattice units) 



^k. ^^K + 0.034(2) - 0.06(2)m^ , (40) 



m^^^ mi + 0.034(2) -0M{2)mi 

™2 



where and ttik are the corresponding pion and kaon masses calculated from domain-wall 
correlation functions. 
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APPENDIX A: MIXED MESON TWO-POINT CORRELATION FUNCTIONS 

Understanding the two-point correlation functions of mixed mesons composed of one 
"Wilsonized" -staggered fermion and one domain- wall fermion, whose interpolating operators 



are given in Eq. (21), requires an understanding of the symmetries of the staggered lattice 
action. This has not been previously worked out before, although pieces of the required 
formalism can be found for example in Refs. 



In the above interpolating field, ^jjdwix) is a domain- wall fermion operator and ^psi^) is a 



'Wilsonized" staggered fermion operator as discussed in the text in Sec. IIIB, and T is one 
of the set of gamma matrices intended to provide the spin and parity quantum numbers of 
a given meson. As we shall demonstrate, because of the peculiarities of the staggered lattice 
action, for all F in 

F = {1, 75 , 7a* , 757a* , If^liy} (Al) 
the corresponding mixed meson interpolating operators have a non-vanishing coupling to 



pions. In particular we will show that Eqs. (25) provides the correct description of the long 



time behavior of the mixed two-point correlation functions. 

We begin with a review of the doubling symmetries of the naive fermion action, which 
are also symmetries of the staggered action and provides a basis for the following discussion. 
Much of the first part of this discussion can be found in Ref. ^85j . 

The free, naive and staggered fermion actions are invariant under the following "doubling" 
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transformations 



S — > S = a^^V(a;) 



V,. 



-7m + ^ 



(A2) 



where in four space-time dimensions g are the 16 elements of the set of ordered hsts with 
up to 4 elements (including the empty set): 

G = {g : g = (/ii, /i2, /^3, At4), /ii < /U2 < /is < /i4} , (A3) 
^5 = n^^ ' with = ^757^ . (A4) 
and the momentum shifts are given by 

(A5) 

0, fi^g 

Here we are using Hermitian Euclidean gamma matrices satisfying the ant i- commutation 
relations, Eq. (18). To realize the doubling symmetry, one introduces 16 momentum space 
spinors. 



[IT. 



g^(fc) = M,^(fc + 7r,), 



(A6) 



which when combined with the periodicity of the integrand over the Brillouin zone, allows 
one to show that the naive action has a 16-fold degeneracy (in 4 space-time dimensions). 



S 



^i-f^-sm{kf,a) + 



m 



d^k 



-./2a W 

d'^k 

-./2a 



tp{k + Tig 



^i7^-sin([A; + 



TTgJ^a) + m 



i){k 



Tin 



^27^-sin(fc^a) + 



(AT) 



To connect the staggered and M^z/sonz^Cfi-staggered fermions with the naive fermions, it 
is useful to use the Kawamoto-Smit transformation of the naive fermion action in position 
space [En]. First one introduces new fields, 

^{x) = n\x)-^{x) , (A8) 
$(x) = ip{x)n{x) , (A9) 
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where 



nix) = iii^,r-/\ (Aio) 



The action expressed in terms of these new fields is then 



S^ = a'Y,Yl '^"(^) k(^) — + H '^o.ix) , (All) 



a 

with 

ry^(x) = (A12) 

In this form, the action is a sum over the four independent spin components of the Dirac 
fermion, $, and one can immediately relate the single component staggered fermions to the 
14^i/sonz2e(i-staggered fermions or the naive fermions by making use of a constant spinor, 

^aix)=C^xix). (A13) 

This also allows us to relate the propagators of the staggered and VFz/sonzzefi-staggered 
fermions which is necessary to understand the mixed meson two-point correlation functions. 
One finds the propagator from x to y, 

D^^{x;y) = niy)I,D^{x;y) n\x) 

= n{y)n\x)D^{x-y). (A14) 

This analysis also allows one to easily derive the behavior of staggered meson two-point 
correlators by converting the various meson sources from the naive fermion form to their 
corresponding staggered form as is found in Ref. [HI]. We now describe the source of time- 
oscillatory behavior. 



1. Time oscillations from doubling symmetry 

For our analysis of the mixed correlators, it will be convenient to examine the fermion 
fields in momentum space with the momentum split up amongst the different corners of the 
Brillouin zone. To do this, we first introduce the subset of G which contains the spatial 
momentum shifts only, Gg- Using the periodicity of the momentum region over the Brillouin 
zone and Eq. (A6), one has 

and similarly for the domain- wall fermion 

The full doubling group is recovered with the product of the spatial doublers and the time- 
doublers, gg and gsgt, with Hg^ = ^(0, 1) and 

M,, =^7574. (A17) 



19 



However, to understand the source of time-oscillations present in the mixed meson corre- 
lators, it is useful to treat the time doublers separately, which we now do. Using the time 
doubling symmetry, gt, one can show that the staggered (and naive) field can be expressed 
as 



— * 



4 ikit 



2tx 



dk 



7r/2a 



2ti 



dk 



■n/2a 



27r 



'Mlq^{k,k,) + 



^irMl^q^^^{k,k,) 



(A18) 



One can immediately see that this field contains both "straight" as well as time-oscillating 
components. Close to the continuum limit, the time-oscillating quark field is of a different 
staggered taste than the straight state. At finite lattice spacing, there are taste-changing 
interactions and so one can not create a source with definite taste. However, we argue that 
at finite lattice spacing, the oscillating states will still arise from this particular taste of 
staggered fermion, or more precisely with a particular change in taste, with the momentum 
shifted to the "time" corner of the Brillouin zone. We will address these issues further when 
we analyze the two-point correlation functions with Dirichlet boundary conditions. 



2. Staggered symmetry and mixed meson correlation functions 

To understand the behavior of the mixed meson two-point correlation functions, we must 
go beyond this, which requires a more detailed analysis of the staggered symmetries. Con- 
structing the staggered fermion action from the naive one amounts to a maximal diago- 
nalization of the doubling symmetry group, G, in which one chooses a maximal subgroup 
H G G, for which the matrices, Mh, defined in Eq. (A4), are all commuting [88] . 

[M,,,Mh>]=0 , y{h,h'}eH. (A19) 

This provides constraining equations amongst the 16 fermion doublers, or "tastes", reducing 
the fermion degrees of freedom, and leading to the relations 

ts{x) = d^'^-^Mhi^s{x) , (A20) 

where 

Mh = ^. (A21) 
To provide specific examples, we will use the set 

H = {gf,:gt, = 0, (3), (1, 2), (1, 2, 3)} , (A22) 



with 

and therefore 



Ao = Ai2 = l , A3 = Al23 = ^, (A23) 
M0 = I4x4 , Ms = ^7573 , Mi2 = -ilil2 , Ml23 = 74 ■ (A24) 
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To understand how the staggered symmetry modifies the analysis, we must return to the 



decomposition of the staggered field in momentum space, Eq. (|A15|). First, we learn that a 

(A25) 



time-momentum space equivalent of the the staggered symmetry, Eq. (A20) is 

^s{k,t) = Mhips{k + 7ih,t). 
Reexamining the staggered field in momentum space, we find 



ijjsix) 



../2a (2vr)3 



h&H gs&Gs/H 



heH gs(^Gs/H 



(A26) 



We then insert this field into the mixed meson interpolating operator, project onto zero 



spatial momentum, and use Eqs. (A15) and (A16) to find 



x 



/•'^/2a d^k - -> + ^ 

E / j:^^i'du,{k,t)T Mlij,{k,t) 
f^^J-./2a (2vr)3 



./2a (2vr)3 



(A27) 



A key observation to make is that the domain-wall fermions do not share the time-doubling 
symmetry that the naive and staggered fermions have. Therefore, the modes in the sum 
over Qs ^ h represent mesons with momentum on the order of the cutoff, p Tt/a and can 
therefore be integrated out of the theory and neglected in the analysis. The general mixed 
meson interpolating field couples to the following sources, 



d^k 



./2a (2vr)3 



iJdw{k,t)T'ips{k,t) + ipdw{k,t) T-/iil)s{k,t) 

+ i^dw{k, t) 17573 ^ps{k, t) + iipdw{k, t) 17271 t) 



(A28) 



To understand the corresponding states this excites, we consider the general mixed meson 
two-point correlation function, 

cl^{t, p = 0) = x) rt x) ^'^^(o, 0) r V'.lo, 0) ) 

X 

= Y.{Dd^{Q, 0; t, x) T D^it, x; 0, 0) ) 

X 

= a' Y^iDdUO, 0; t, x) r D^{t, x; 0, 0) n\t, x) ) 

X 

= a^Yl ^(^) < ^'^-(0' 0; ^x(0' 0; x) r n\t, x) ) . (A29) 
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with 

e(x) = (-l)(^i+^2+x3+t)/a^ (A30) 

and ( ) denoting a sum over the color and spin indices and an average over gauge configu- 
rations. Strictly speaking, this equation, and similar ones that follow, are only true in the 
limit of an infinite number of configurations. The spin factor, F x)r''" for the various 
r will determine which staggered phases and quantum numbers a given source excites, for 
which the following relations are useful 

75 fit(t,x) 75 = e(a;) ^^^(t,x) 

7^ n^{t, x) 7^ = 7]^ix) C^ix) n\t, x) (A31) 

with the remaining staggered phase factors given by 

Viix) = 1, V2ix) = i-ir\ vsix) = r],{x) = (-l)-^+-2+-3, 

Ci(a;) = C,{x) = C3(x) = (-l)-^ Ux) = l. (A32) 



In addition to determining the states excited by the sources in Eq. (A28), we can also verify 
that the time-oscillating partners related by a change in the staggered quark taste, from 
Eq. (A18),^*^ otherwise have the same quantum numbers, for which the Dirac structure is 



modified by 7475, Eq. ( |Al7l ). We find 



CJs (t) = 5^ ( D,^ (0, 0; t, x) (0, 0; t, x) fi^t , x) ) , (A33a) 

X 

Qtit) = i-iy+'a^ D,^(0, 0; t, x) D*^iO, 0; t, x) Q\t, x) ) , (A33b) 

X 

CJf ^(t) = a'J2 Viix){DdUO, 0; t, x) (0, 0; t, x) n\t, x) ) , (A34a) 

X 

Cl{t) = i-iy+'a' J2 ^4 W ( ^^-(0, 0; t, x) ^^(O, 0; t, x) n\t, x) ) , (A34b) 

X 

Qm = a' Y,(-^T' ( ^'^-(O' 0; ^' ^) ^x(0' 0; x) n\t, x) ) , (A35a) 

X 

CJf^W = (-l^'a3^(-l^+^■^r74(x)(D,^(0,O;t,x)Z};(0,O;t,x)^]t(^,x)), (A35b) 

X 

Cjr(t) = Y.^-irv^ix) ( D,^{0, 0; t, x) ^^(0, 0; t, x) n\t, x) ) , (A36a) 

X 

CJrW = (-l)*+V^(-irr/4(a;)(Dd^(0,0;t,x)D*(0,0;t,x)(^t(i^x)), (A36b) 

X 

explicitly verifying both the time-oscillatory nature of the correlation functions related by a 
time-doubling transformation. 

This can be equivalently understood working in position space and applying the time-doubhng transfor- 
mation, Ipsix) ipsix) and determining the change in the two-point correlation functions. 
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P. = '/'dwYj'/'s 
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(a) 



(b) 



FIG. 2: Correlation functions from mixed meson sources. In Figure|2]^o), we plot the correlator for 
the pseudo scalar source, iTyg = '4'dwl5ips- In Figure [2|&) we plot the mixed "rho" meson correlator, 
Pz = V'dwiTsV's- The long time behavior of both correlators is dominated by single pion states, as 
demonstrated by the best fits to both with Eq. (32). See Sec. |A] for details and Table IV for results. 



3. Pions from p's, ao's and other meson sources 



We can now show that pions couple to all of the naive mixed meson sources with Dirac 



structure given by Eq. (22), by inserting T = {1,75,73,7573,7172} into Eq. (A28), and 
using Eq. (A18), one can verify that all of the TT^g^x) interpolating fields couple to all the 



sources given in Eqs. (A33a)-(A36b). Because of the mixing of spin and taste degrees 



of freedom in the staggered action, computing the two-point correlation function given in 



Eq. (A33a), one is not only exciting states with pseudo-scalar quantum numbers. Using 



the symmetries of the staggered action, we have shown that the source for this correlator 



couples not just to Eq. (A33a), but that Eq. (A33a) is just one component of a larger 



matrix of correlation functions excited by the naive pion source, for which Eqs. (A33a)- 



(A36b) are diagonal entries. The naive pseudo-scalar source will have an overlap with the 



propagating pion, as one term of this source has the quantum numbers of an odd parity 



scalar, but likewise from Eq. (A28), so will all the naive mixed meson sources. Therefore 



we have shown that Eq. (25) provides the correct description of the long-time behavior of 



all the mixed meson interpolating fields created with (anti) periodic boundary conditions 
considered in this work (anti periodic boundary conditions for the fermions and periodic for 
the mesons). In principle, one can determine what linear combination of sources, vr^^ will 
have an overlap with only the pion and other mesons by studying the matrix of correlation 
functions, disentangling the eigen-states excited by these sources. This however, is beyond 
the scope of this work. 



In Table IV, we present the results of fits to the long time behavior of the various mixed 
meson sources with the full volume mOlO propagators and compare them to the results of 
fitting the naive mixed pion source. In Figure [2| we display the best fit plots and results 
for the two-point correlation functions created with the "pion" and "rho" interpolating 
operators. In Figure |3] we plot the best fit and results for the correlation function created 
with the F = 74 interpolating field, which demonstrates an even more dramatic example 
of the time-oscillatory nature of the mixed meson correlators. In this case, the oscillating 
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FIG. 3: Correlation function created from a mixed Oq^ source. In Figure [3][|a) we see the dramatic 
oscillatory effects which can occur when the oscillating amplitude is larger than the straight am- 
plitude for degenerate or nearly degenerate masses. In Figure |3|6) we show two different fits to it 
described in App. |A 3 



TABLE IV: Mixed meson from p and aO and b mesons. We collect our results of fits to the long 
time behavior of the various two-point correlation functions created with the listed sources and 
compare the to that of the mixed pseudo scalar source. 



Ensemble 


source 




fit range 


/ d.o.f. 


Q 


mOlO(FV) 


'>Pdw 75 Tps 


0.2851(31) 


6-31 


20 / 23 


0.65 


mOlO(FV) 




0.283(19) 


18-31 


12 / 11 


0.33 


mOlO(FV) 




0.289(04) 


8-31 


17.6 / 21 


0.68 


mOlO(FV) 


V'dto 73 V's 


0.284(08) 


9-31 


17 / 20 


0.68 


mOlO(FV) 


V'dto 72 V's 


0.285(09) 


11-29 


8.9 / 16 


0.92 


mOlO(FV) 


i'dw 7172 V's 


0.271(18) 


17-32 


8.0 / 13 


0.84 


mOlO(FV) 


i'dw 7375 V's 


0.282(23) 


21-27 


2.4 / 4 


0.67 



exponential has an amplitude roughly two orders of magnitude larger than the straight 
exponential as seen in Figure |3](a). In Figure [3|6) we have actually plotted two different fits 
to the correlator, one with a forced degeneracy between the straight and oscillating state 
and the second letting the two mass parameters be independent, 



Ci{t) 



A + Bcosinrit) 



-mt 



C2{t) = Ae""' + Bcos{Tmt)e 



-Mt 



(A37) 
(A38) 



and we find consistent masses in both cases. In Table IV, we collect the results of our deter- 



mination of the mixed pion mass from all of the two-point correlation functions presented 
in this Appendix. 
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4. Mixed mesons and Dirichlet Boundary Conditions 



In this appendix, we demonstrate the equivalence between the time shift symmetry and 
the taste-changing transformation which gives rise to the time-oscillatory states. This makes 
it plausible that the Dirichlet boundary conditions, which break the time-shift symmetry, 
give rise to the observed splitting of the straight and oscillating states observed in the mixed 



meson correlation functions, Eq. (27). 

In the continuum limit, one can associate the time oscillating field with a specific taste 
of staggered fermion field, most readily seen in momentum space. Furthermore, using the 
momentum space time axis-reversal transformation for the staggered fermion fields [901 E] , 
one can show that under time-reversal, the staggered field undergoes an identical taste change 
as that from the time-doubling symmetry. With time-anti periodic boundary conditions for 
the fermions, both of these transformations are symmetries of the staggered action, and 
therefore we conclude that the straight and oscillating states in the two-point correlation 
functions must be degenerate. 

To show that the transformations of time-doubling and time axis-reversal correspond to 
the same taste change of the staggered field, it is simplest to work with the momentum space 
fermions. To do this we construct a 16 component staggered field with components defined 
by the fermions at the different corners of the Brillouin zone, 

Qg{k) = x{k + -K,), (A39) 



with Tig defined in Eq. (A5). Then using the symmetries of the staggered action, one can 



show that under a time axis-flip, this fermion transforms as 

g(A;) ^r4r5S4S5Q(fc), (A4o) 

where the set of 16 x 16 matrices, F and S can be found for example in Refs. [901 E] • If we 
begin with the staggered field in the central corner of the Brillouin zone, Qo{k), then one 
can show 

T^T,E,E,Qo{k)=Q^{k), (A41) 



where tt^ = - = (0, 0, 0, 1). From Eq. (A18), we recognize that this is the exact taste change 



which occurs with the time-doubling symmetry. Furthermore, if we apply this time axis-flip 



to all of the components of the staggered field given in Eq. (A26) which are related by 



the constraining equation, Eq. (A20), then one can show that the taste change from the 
time-doubling transformation and the time axis-flip transformation are identical. 

This makes it plausible that the use of Dirichlet boundary conditions, which break this 
symmetry, give rise to the observed mass splitting between the straight and oscillating states 
observed in the mixed meson two point correlation functions. We reiterate that this is not an 
issue for the mixed action calculations to date [IlEllllElEl[ZllHl[ni[IDllIIllI2llI5llIS],as they 
only impose Dirichlet boundary conditions on the valence domain-wall fermions. Our use of 
the Dirichlet boundary conditions for the staggered propagators used in the construction of 
the mixed meson two-point functions was one of convenience. 

In Table |V| we present our analysis of the quantity 

(aMr - {arhl'f , (A42) 
which shows no discernible dependence upon the pion mass. 



25 



TABLE V: Mixed meson and Dirichlet boundary conditions. We present the mass splitting of the 
two mixed pions arising in a MA scheme with DBCs. This mass spHtting shows no dependence on 

the ("luark masses within our |)rocision. 



Ensemble 


quantity 


result 


fit range 


/ d.o.f. 


Q 


m050(416 DEC) 


TT 


0.532(3) 

V / 


6-13 


4.95 / 4 

/ 


0.292 


m050(416 DEC) 




0.500(2) 


6-13 


4.95 / 4 


0.292 


m040(349 DEC) 


aM^' 


0.489(5) 


6-13 


6.39 / 4 


0.172 


m040(349 DEC) 


am^^ 


0.457(2) 


6-13 


6.39 / 4 


0.172 


m030(564 DEC) 




0.443(4) 


6-12 


9.48 / 3 


0.02 


m030(564 DEC) 




0.409(2) 


6-12 


9.48 / 3 


0.02 


m020(222 DEC) 


aM^' 


0.393(6) 


6-12 


1.64 / 3 


0.650 


m020(222 DEC) 




0.352(3) 


6-12 


1.64 / 3 


0.650 


m010(447 DEC) 


aM^' 


0.324(11) 


6-12 


4.51 / 3 


0.21 


m010(447 DEC) 




0.282(3) 


6-12 


4.51 / 3 


0.21 


Ensemble 


mOlO 


m020 


m030 


m040 


m050 




0.025(08) 


0.034(11) 


0.029(03) 


0.031(04) 


0.033(04) 
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